Evidence of growing spatial correlations during the aging of glassy glycerol 
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We have measured, as a function of the age t a , the aging of the nonlinear dielectric susceptibility \3 
of glycerol below the glass transition. Whereas the linear susceptibility can be accurately accounted 
for in terms of an age dependent relaxation time T a (t a ), this scaling breaks down for X3> suggesting 
an increase of the amplitude of X3- This is a strong indication that the number N corr of molecules 
involved in relaxation events increases with t a . For T — 0.96 x T g , we find that N corr increases by 
~ 10% when t a varies from Iks to lOOks. This sheds new light on the relation between length scales 
and time scales in glasses. 



An old plastic ruler under tension has a longer length 
than a newly made one [l| . This is a striking illustration 
of the aging phenomenon, the hallmark of the physics 
of glasses. The physical properties of an aging system 
depend on the time t a elapsed since the material has 
fallen out of equilibrium; i.e., since the glass transition 
has been crossed. Understanding aging is of paramount 
importance [l|, Q , both from a fundamental and a prac- 
tical point of view (many daily-life materials do in fact 
age). Yet, there is no universally accepted theoretical 
description of the basic mechanisms of aging, although 
many scenarii have been proposed (see [2H6|])- 

One of the most distinctive features of aging is the in- 
crease of the physical relaxation time r Q with the age t a 
[2> EL 04l3- I n t ne case of spin-glasses, this increase has 
been rather convincingly attributed to the growth of the 
number Nc 0rr of cooperatively relaxing spins. Both simu- 
lations [l3( and experiments [l4( are compatible with this 
scenario, and allow one to estimate the dynamical growth 
law N corr (t a ). The situation is much less clear for most 
other gl assy systems, either experimentally or numeri- 
call v I13l . ll5| . In fact, while there are only two simulations 

, reporting the growth of a dynamical correlation 
length during the aging of model glasses, and few ana- 



lytical studies [13, ll 



j, there is to our knowledge no 
available experimental result for real glass-formers. The 
last decade has witnessed an outburst of activity on dy- 
namical heterogeneities and on the determination of the 
size N corr of dynamically correlated molecules in glasses 
0, HI, but almost all these studies have been confined 
to equilibrated systems. Whereas a compelling positive 
correlation between N corr , eq (T) and the equilibrium re- 
laxation time r a (T) has been established above the glass 
temperature T g (see [H and refs. therein), its aging coun- 
terpart has not been investigated experimentally. The 
aim of the present study is to extend to the aging regime 
the experimental determination of N corr that relies on 
the cubic nonlinear dielectric susceptibility [2(J 21 1. We 



As argued in 19], non-linear susceptibilities are the 
ideal gambits that elicit the growth of amorphous order 
in glassy systems. Whereas linear susceptibilities (dielec- 
tric, magnetic, elastic, etc.) are blind to amorphous order 
and dynamical correlations, the equilibrium cubic nonlin- 
ear dielectric susceptibility X3 01 deeply supercooled glass 
formers is given, at temperature T, by [1 91 ] : 



X3(w,T) « Z(T)N t 



corr.eq 



(T)H(uT a (T)) (1) 
(2) 



will report, for the first time, clear experimental evidence 
of the growth of the size of dynamically correlated regions 
during the aging of glycerol - a prototypical glass former. 



where uj = 2nf is the angular frequency, fcs the Boltz- 
mann constant, eo the vacuum permittivity, a 3 (T) the 
molecular volume, and A\i(T) = Xi(w = 0, T) — xi(cj 3> 
r^ 1 ) is the contribution to the static linear susceptibil- 
ity of the degrees of freedom associated with the glass 
transition. In Eq. ((T|), H(u), with u = ujT a {T) = 
f/fa(T), is a complex scaling function which goes to 
zero both for small and large arguments, and peaks 
in-between. Eq. ([1]) can be fully justified within the 
Mode-Coupling Theory of glasses [221 ]; it has been con- 
firmed experimentally in [2(J [2l|, and used to extract 
precise estimates of N corr ^ eq (T) in equilibrium. In the 
aging regime, it is natural to conjecture [l|| that the 
above expression remains valid with r a (T) — > T a (t a ) 
and N corr>eq (T) — > N corr (t a ), therefore allowing one to 
infer information about the growing of N corr during ag- 
ing. Strictly speaking, such a simple substitution is too 
naive: one expects on general ground that (a) the scaling 
function % should also be replaced by a different scal- 
ing function H; and (b) the prefactor Z(T) might itself 
acquire an age dependence: the value of both Axi and 
a could evolve with age, and the temperature T should 
in principle [l|| be replaced by an "effective" temper- 
ature T e ff(t a ) that encodes the possible deviations to 
the equilibrium fluctuation-dissipation theorem (23l - [25| . 
However, our experiments are "weakly" out of equilib- 
rium, since they reach equilibrium eventually. In this case 
we expect that the scaling assumption Eq. (fT]) general- 
ized to the aging regime, with H = H and Z{t a ) = Z(T), 
holds to a very good approximation. Our strategy will 
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For comparison, we also plot the equilibrium data at 
T„ + 6K obtained in [20|, |2l| which shows 
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FIG. 1: (Color online) Equilibrium values of the modulus and 
the phase of X3, measured in glycerol for three temperatures 



around T„ 



188K. 



therefore be the following: (i) since the linear suscep- 
tibility does not depend on N corr , its age dependence 
should only come from that of T a (t a ). Indeed, we will 
establish that t a ) ~ G"(wT a (t a )), where Q" is the 

equilibrium scaling function (see Fig. [3]). This allows us 
to determine T a (t a ) directly; (it) by waiting long enough 
(i.e. t a — 200ks) we measure the equilibrium non-linear 
susceptibility X3(oj,T) for various frequencies (see Fig. 
[3]), thereby allowing one to obtain the scaling function 
T-L(u); (Hi) inserting these informations into Eq. ([ij in 
the aging regime, we can deduce the age dependence of 
Ncorr (up to the assumption that Z(t a ) = Z(T), see Fig. 

ED [2g. 

Experiments. Ultrapure glycerol was purchased from 
VWR and placed in our dielectric setup described in Refs 
[20I 21 , 27 1 . Glycerol was the dielectric layer of a ca- 
pacitor made with stainless steel electrodes separated by 
a 8.25^m thick Mylar© ring. All the aging quantities 
were measured with the same T quench: The sample 
was first set to 196K ss T g + 8K (where r a « Is) during 
lhr, then it was cooled, without any undershoot, to the 
working temperature T = 180. IK (or T = 182. 7K) in 
1.8ks, and finally T was kept constant within a ±70mK 
interval during 200ks. We used a high harmonic purity 
a.c. field of amplitude < 3MVrms/m to measure, sepa- 
rately, X3(w,i ) as well as Xi(w,t ) and Xi(3w,t ) " see 
[2^ . Once X3( w ?0 is known, the key quantity is, ac- 
cording to Eq. (|TJ|, X 3 (uj,t a ) defined as: 

X 3 (u,t a ) = X3 ^;\ a) « N corr (t a )U(LOT a (t a )) (3) 

The equilibrium values of Xj, yeq (u>,T), obtained after 
the end of aging, are plotted as a function of f/f a . eq 
in Fig. CD for T = 180. IK and for T = 182.7K (f a 
is denned as the peak frequency of tfl eq (w,T), |29{). 



that the qualitative trend already found above T g [20, 21] 
holds also below T g , i.e. for fixed / / f a ,eq(T) in- 

creases when T decreases, [30j. 

Scaling analysis. During aging both |xi| and IX3I de- 
crease for a given cj, because the dielectric spectrum 
shifts to lower frequencies @-[l2j- This is illustrated by 
the series of symbols on Fig. [2] For not too deep quenches 
0, such as ours, and for our limited range of 3.5 decades 
in frequencies, the shape of the spectrum of xi' is not 
expected to change during aging except for an overall 
scaling factor. This is confirmed by Fig. [2] where the 
Xi (to, t a ) data for all frequencies and all ages can be very 
accurately reproduced by the equilibrium susceptibility 
Xi eq (w,T = 180. IK), up to a rescaling of the frequency 
by a factor x(t a ) = f a (t a )/fa,eq (see the dotted lines in 
Fig. [2]obtained by adjusting the factor x(t a ) for each t a ). 
We have also checked that our values of x(t a ) are close 
to what is predicted by the ansatz introduced in Ref. [§] , 
[3lj ]. We have checked that the very same x(t a ) factor 
also allows us to rescale the x'i( w j ta) data onto the equi- 
librium curve. Note that since the x" (<*>, t a ) are not pure 
power-laws in frequency, horizontal and vertical shifts (in 
log-log) are not equivalent. Hence, the accurate rescaling 
of Fig. [2] implies that the amplitude of Xi ( w ; t a ) does not 
depend on the age. A finer look at the rescaling suggests 
that this amplitude is constant within a 1% uncertainty 
range, and if anything, decreases with age. This is impor- 
tant for the discussion of the Z(t a ) factor in Eq. (p} that 
includes Axi{t a ) to which Xi( w )0 is proportional. To 
estimate the difference between Z(t a ) and Z(T) we in- 
voke the "Active" temperature Tf ict (t a ) 0, @, 0, H3| (not 
to be confused with the effective temperature T e //) de- 
fined such that Xi (<*>:*<«) = Xi, eg (^> Tf ic t(t a )). This phe- 
nomenological recipe leads to Tf ict (t a > 0.3ks) < T+2K. 
By extrapolating the T dependence of Axi above T g , we 
estimate that Axi(t a ) may increase by at most 1.1% dur- 
ing aging. The order of magnitude is similar to the one 
suggested by the rescaling analysis above, albeit with an 
opposite sign. Similarly, we estimate that a 3 might de- 
crease by ~ 0.2% during aging. Finally, close to 180K, 
T/T e ff(t a ) was found to increase in glycerol during ag- 



ing, by ~ 2% according to [24J, but by at most 0.7% 
according to the recent work of Ref. (2||. Altogether, we 
conclude that Z(t a )/Z(T) remains very close to unity, 
with a probably much overestimated maximum increase 
of 4% during aging. This effect is therefore smaller than 
the ~ 12% increase of N corr (t a ) that we infer from our 
analysis below. 

Our central experimental result is summarized in Fig. 
[3] where we now show both \\ ( w > ta) and \X 3 (uj,t a )\ as a 
function of f / f a (t a ) for Iks < t a < 200ks. As expected 
from the results of Fig. [21 xi'^^a) collapses very well 
onto the equilibrium curve. However, this collapse is not 
observed for \X3(ui,t a )\ (Fig. [31 filled triangles and left 
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FIG. 2: (Color online) Aging of the out-of-phase susceptibil- 
ity Xi(w,*o) of glycerol at T = 180. IK (filled symbols), in 
log-log coordinates. The thick (resp. thin) solid line corre- 
spond to the equilibrium spectrum at T = 180. IK (resp. at 
T = 182. 7K). The dotted lines superimposed to the filled 
symbols are obtained by translating horizontally the 180. IK 
equilibrium spectrum by a factor x(t a ) (see text). Inset: Re- 
sulting curve for the scaling factor x(t a ) = fa(t a )/ fa,eq- 



FIG. 3: (Color online) T = 180. IK. Solid curves: equilib- 
rium quantities vs f/f a ,eq- Small symbols: aging quanti- 
ties vs f/f a (t a ) -the value of / labels each data set-. Note 
the collapse of the aging and equilibrium curves for x" (right 
axis). This contrasts with what is observed for \Xs\ (left axis) 
where the aging values are systematically below the equilib- 
rium ones. This reflects the increase of N corr (ta) during ag- 
ing; see text, Eqs. (|3I4|I and Fig. 0] 



axis). The rightmost points of these series of triangle cor- 
respond to the equilibrium values \X3 teq (w,T)\ and are 
singled out as large black squares. The thick line joining 
these black squares is an interpolation that corresponds 
to the equilibrium value of \Xs te q(u},T)\ for intermediate 
frequencies. At variance with the good superposition ob- 
tained for x"i Fig- El reveals that, for a given f/f a (t a ), 
the value of \X3(uj, t a )\ is systematically below the cor- 
responding value at equilibrium. This is exactly what is 
expected from Eq. (J3J: the ratio between these two val- 
ues should be equal to N corr (t a ) / N corr ^ eq (T) , and should 
thus increase with age, precisely as observed in Fig. [3] 
Defining the vertical logarithmic distance <5, in Fig. |3l 
as: 



S(u,t a ) 



Z(t a )X 3 (uf a (t a ),t a ) Z 

(^a)-^ corr (^a) 



Z(T)X 3 , eq (uf a , eq ) Z(T)N corrte q(T) ' 



(4) 

we obtain the last equality if Eq. ([3]) holds, in which 
case S(u, t a ) should be independent of frequency. With 
Z(t a ) / Z(T) « 1 justified above, we conclude that 6(u,t a ) 
directly measures N corr (t a )/N corr ^ q (T). 

The values of 5(u,t a ) are plotted in Fig. H |33] . We 
indeed observe that, up to the precision of our measure- 
ments, 5(u,t a ) does not depend on frequency. This is 
an important consistency test of our scaling assumption, 
Eqs. (11131) . From Fig. |H we deduce that S(u,t a ) increases 
by ps 12% when t a increases from Iks to lOOks. Since the 
ratio Z(t a ) / Z(T) increases by at most 4%, we interpret 
the data of Fig. 2] as giving the first experimental evi- 
dence that the size of the dynamically correlated clusters 



increases with the age in a glass former, see [34|. The in- 
crease of N corr (t a ) during aging can be approximately 
accounted for by extending the observation made in 



the temperature dependence of N, 



corr,eq 



deduced 

from non-linear susceptibility measurements can alterna- 
tively be obtained as: dN correq /dT w l.bd(Txr) I dT , 
where Tyt = Tx max u \ [d(xi,e q (u,T)/A X i)/dT] |, see 
21, 35 371. We now surmise that this can be extended to 



the out-of-equilibrium regime by simply translating the 
fait g) dependence of Fig. [2] in terms of Tfi ct (t a ) (see 
[3q). This heuristic procedure leads to the solid curve 
shown in Fig. |H which is indeed close to the values of 
5{u,t a ) directly drawn from our experiments. This sug- 
gests that it mig ht be possible to extend the theoretical 
work of [35I l36| to aging, and get a simplified way of 
estimating N corr (t a ) using linear susceptibilities. 

Time and length scales. Finally, we take advantage 
of the wide range of time scales over which the evolu- 
tion of N corr has been measured to revisit one of the 
most crucial aspect of glassy dynamics, namely the re- 
lation between time and length scales. Within the Ran- 
dom First Order Transition theory 3{| 4(|, one expects 
ln(r Q /ro) = T^p S /fcsT, where To is a microscopic time 
scale, T a typical molecular energy barrier, ip$ the point- 
to-set correlation length 0, [39|, which sets the size of 
the clusters that must rearrange cooperatively for the 
system to relax, and ip the so-called barrier exponent. 
In Wolynes' version of RFOT, T = kIcbT where k is a 
number that depends weakly on molecular details, and 
-0 = 3/2 [40]. In order to compare with our results, 
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FIG. 4: (Color online) T = 180. IK. Values of S(u,t a ) = 
Z(t a )N corr (t a )/Z(T)N corr (eq) « N corr (t a )/N corr (eq) ex- 
tracted from Fig. [3] S(u,t a ) is found to be independent 
of frequency, as predicted by Eqs. 1)314(1 . The solid line 
is the estimation based on phenomenological Active tem- 
perature Tfict{t a ), see text. Inset: Estimation of the ex- 
ponent ip (see text) using all aging and equilibrium data, 
starting from 204K where T a ,eq — 20ms. The :r-axis is 
R(tc) = iV corr (r C( )/JV corT . >eg (20ms) with r a = r a (t a ) or r a ^ q . 
For the y-axis, we chose either T = To which amounts to 
V = ■u(T Q ,T)/i>(20ms,204K) and yields ip « 1; or T = nk B T 
which amounts to W = tu(T a )/uj(20ms) and yields ip ~ 3/2. 



one should postulate that the size of dynamically cor- 
related clusters N corr is proportionnal to Cpg. This re- 
lation is not unreasonable, but sharp theoretical argu- 
ments are still lacking to relate unambiguously "coop- 
erative" regions to "dynamically correlated regions" . In 
any case, we collect all our past and present data in the 
inset of Fig. |H where we plot v(r a ,T) — Tln(r Q /r ) 
and w(r a ) — ln(r a /ro) as a function of N corr , where 
T a = T a (t a ) or r a , eq . In Fig. V and W correspond 
respectively to T = To independent of temperature and 
to T = nksT. We fix the value of To to 10 _14 s, leav- 
ing ip as the only free parameter. To our surprise, we 
find that the best choice for ip within the first hypothesis 
(T = To) is ip ~ 1, which is the value found numerically 
in [4l|, whereas in the second hypothesis (T = k/cbT), 
we find ip ~ 3/2, as predicted by Wolynes et all Our 
data is compatible with both hypotheses j42j, although 
slightly favoring the first one, in particular in the aging 
regime (see Fig. [4}. Note that a factor 10 on To changes 
the value of ip by ~ 10%. 

Conclusion. We have reported the first direct obser- 
vation of the increase of N corr in the aging regime of a 
structural glass. For glycerol at T = 0.96T g , our quench 
protocol yields an increase of N corr by ~ 10%, which 
lasts 



glasses. Our study opens a new path for studying aging 
in many other systems. It could also be extended to the 
more complicated thermal histories designed to probe the 
memory and rejuvenation effects @, @|- Monitoring the 
behaviour of N corr in these experiments should shed a 
new light on these phenomena. 
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lOOks 34]. These results deepen our micro- 



scopic understanding of aging and give precious informa- 
tion about the relation between time and length scales in 
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SUPPLEMENTARY INFORMATION: 

For the sake of completness, we give here the main in- 
gredients of the 3w measurements in the aging regime. 

(3) 

Note that the quantity defined as X3 i- n Eq. (6) below 
has been noted xs in the main letter, to simplify the no- 
tations. 
As explained in Refs. 



applied onto a dielectric liquid, the macroscopic polari- 
sation P can be expressed as : 



X i(t-t')E(t')dt' (5) 

X3(t-*i,t-*2>*-t3)-E(ti)x 

) 

xEiQEit'^dt^dt'zdt'z + 



where the function Xi(t) corresponds to the experi- 
mental macroscopic linear response while X3(*ij^2,^3) is 
the experimental macroscopic nonlinear response. This 
expression is valid as long as the non linear terms are 
small, i.e. as long as, for any integer k, \x2k+iE 2k+1 \ <C 
\X2k-iE 2k ~ x \. This is why Eq. (5) is restricted to the 
cubic response, and neg lects higher order terms. 

It is shown in ref. [27], that for a field E(t) = Ecos(wt) 
one gets, from Eq. (5): 



E 



\Xi\ cos(cjt - 61) 
+3/4£ 3 | X < j 1) |cosM-4 1) ) 
+ l/4£ 3 | X f|cos(3^-4 3) )- 



(6) 



The time dependent polarisation amounts to an elec- 
trical current given by 



I(t) = S 



8P(t) 
dt 



(7) 



where S is the surface of the electrodes. Inserting Eq. (6) 
into Eq. (7), one finds that I(t) is the sum of a current 
I(lu),t), oscillating at the fundamental frequency, and of 
a current J(3w, t) oscillating at 3w. As for the field range 
E < 3MVrms/m involved in our experiments the con- 
dition mentionned above \x2k+iE 2k+1 \ <. \x2k-iE 2k ^ 1 \ 
is satisfied, one gets firstly that |/(3w,t)| <C |/(lu;,t)|; 
secondly that the value of I(lui, t) is fully dominated by 
Xi and thus I(lu),t) can be, to a very good approxima- 
tion, analysed by using the usual framework of complex 
admittance y(uj). The small third harmonics current is 
given by: 



I(3u,t) = ^e SoJxi 3) (io)E 3 cos(3o;i + | - A., 



(3) 



) (8) 



21|, [27(, when a field E(t) is 



This quantity is so small that carefully designed elec- 
tronic setups must be used to avoid to mix the sought 
I(3lu, t) with the nonlinear imperfections of the voltage 
source and of the voltage amplifier 27J] . It was shown in 
Ref [27J that two kinds of setups can be used: either a 
"two samples bridge" involving two samples of different 
thicknesses -see the inset of Fig. [S}; or a "twin-T notch 
filter" , see the main part of Fig. [S] 
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The two samples bridge is a technique measuring a dif- 
ferential voltage V rn = Vthin - Vthick and relying on the 
"balancing relation" ensuring that V m (lui) = 0: this hap- 
pens provided one has z thlck y thlck = z thm y thm where y 
is the admittance of one sample and z is the impedance 
relating the sample to the ground. An important fea- 
ture of the two samples bridge is that once the balancing 
condition is met at Iuj it is also met at any other fre- 
quency. Thus the balancing condition enables, at the 
same time, to suppress the contribution coming from 
the nonlinear character of the input voltage amplifier - 
since VthinO-w) - Vthick (lu) = 0-; and to suppress the 
3w spurious component of the source -since the linear 
response of the samples cancels at any frequency-. The 
two samples bridge is thus the most efficient technique 
to get J(3cj,i), provided one has a way to check that 
the condition Vthin(lw) — VthickO-w) = remains true 
during the I(3u), t) acquisition. This condition might in- 
deed not remain true in the case where some uncontrolled 
slight disymetry between the two samples happens, such 
as the one resulting from a slight difference in the tem- 
perature of the two samples. Fortunately, when mea- 
suring at equilibrium, one varies the field E: one thus 
checks all along the 3w acquisitions that /(3a;, t) is ac- 
curately -up to ±1%- proportional to E 3 , which ensures 
that the condition Vthin{lw) — Vthj C fc(lw) = is met - 
enough- during all the acquisitions. However, in the ag- 
ing case, one monitors I(3w, t) as a function of the age t a , 
for a constant E. Thus, if some violation of the condition 
V t hin(lu) - Vthick(lu) = happened, it might pollute the 
age dependence of the 3w response of the samples, and 
we would have no way to correct this imperfection. 

This is why we have decided to work with only one 
sample and to use the "Twin-T notch filter" -see Fig. [5J-: 
its transmission coefficient 7/tw(w) is smaller than 10~ 4 
for the frequency /o = 1/(2ttRC) and of order 1 for 3/o- 
Therefore, by choosing the components R, C so as to set 
2tt/o = lw, we are absolutely sure that the lw voltage at 
the input of the Lock-in amplifier is small enough dur- 
ing the 3w acquisitions in the aging regime. Besides, by 
setting the voltage source to one of the few values where 
the DS360 voltage source is nearly perfectly harmonic, 
one gets a setup where the 3w spurious contribution of 
the source is nearly negligible -this very little spurious 
contribution can, of course, be easily measured and sub- 
tracted from the measured signal-. This is why we have 
choosen the twin T notch filter for our 3w measurements 
in the aging regime. 

At 3lo the sample is equivalent to a pure current source 
I(3ui,t) with an impedance Z{Zw) = l/[y(3w)] placed in 
parallel. Neglecting for simplicity any remaining spurious 
contribution of the DS360 source, one gets for the voltage 
V meas {3ui, i) measured by the lock-in amplifier: 

V meas (3Lj, t) = a(3uj)Z{3uj)I{3uj, t) (9) 




FIG. 5: Main figure: Electronic setup used for the measure- 
ments of x^i^yta) in the aging regime. Note that the quan- 
tity defined as X3 3 ' in Eq. (6) of this Supplementary Material 
has been noted X3 i n the main letter, to simplify the nota- 
tions. The "twin-T" notch-filter damps the response at lu by 
a factor larger than 10 4 . Inset: Two samples bridge, see ref. 



where a(3w) is the global transmission coefficient be- 
tween the sample and the Lock-in. In the simplest case 
where (r, \Z\) <C (R, l/[C3w]) the coefficients multiply 
and one gets a(3w) = Tfutfiu) X r/[r + Z(3v)]. In any 
case, a(3w) is directly measured by setting the funda- 
mental angular frequency of the source to f2 = 3w, and 
by using a(fi) = V meas (Vl)/V s 

o n rce 

(fi). 

Eq. (9) is written at equilibrium, when all the in- 
volved quantities no longer depend on the age t a . At 
equilibrium, we have of course carefully checked, both 
at T = 180. IK and at T = 182. 7K, that V mea3 (3uj,t) 
is proportionnal to the cube of the voltage source V 3 . 
In the aging case, all the susceptibilities of the sample 
-linear and nonlinear- depend on the age t a . As a result, 
all the quantities involved in Eq. (9) depend on the age 
t a . This is why, by repeating for each quantity the very 
same T quench, we have separately measured the age de- 
pendence of I(3uj,t,t a ) but also of Z(3u),t a ), Z(lw,t a ), 
a(3w,t ), and of V ap pi(luj,t,t a ) = V s (lu),t)-V a (luj,t,t a ) 
where V app i(lu!, t, t a ) is the voltage applied onto the sam- 
ple at age t a and at time t -see Fig. [5J-. We have made 
all the possible consistency checks -for example, given 
R, C, Z(lu>, t, t a ) one can predict the age dependence of 
V app i(luj, t, t a )-. A single T quench lasts 200ks at 180. IK, 
and for each of the 5 different frequencies ranging from 
12mHz to 11Hz, we have measured the age dependence of 
the 5 different quantities mentionned above. As a result, 
with all the cross-checks, the data acquisition, for 180. IK 
and 182. 7K altogether, took a bit more than 10 4 ks, i.e. 
a bit more than 4 months. 



